Reply to Barton et al: signatures of natural selection during the Black Death.

Barton et al.1 raise several statistical concerns regarding our original analyses2 that highlight the challenge of inferring natural selection using ancient genomic data. We show here that these concerns have limited impact on our original conclusions. Specifically, we recover the same signature of enrichment for high FST values at the immune loci relative to putatively neutral sites after switching the allele frequency estimation method to a maximum likelihood approach, filtering to only consider known human variants, and down-sampling our data to the same mean coverage across sites. Furthermore, using permutations, we show that the rs2549794 variant near ERAP2 continues to emerge as the strongest candidate for selection (p = 1.2×10-5), falling below the Bonferroni-corrected significance threshold recommended by Barton et al. Importantly, the evidence for selection on ERAP2 is further supported by functional data demonstrating the impact of the ERAP2 genotype on the immune response to Y. pestis and by epidemiological data from an independent group showing that the putatively selected allele during the Black Death protects against severe respiratory infection in contemporary populations.

Helical coiling of metaphase chromatids.

Chromatids of mitotic chromosomes were suggested to coil into a helix in early cytological studies and this assumption was recently supported by chromosome conformation capture (3C) sequencing. Still, direct differential visualization of a condensed chromatin fibre confirming the helical model was lacking. Here, we combined Hi-C analysis of purified metaphase chromosomes, biopolymer modelling and spatial structured illumination microscopy of large fluorescently labeled chromosome segments to reveal the chromonema - a helically-wound, 400 nm thick chromatin thread forming barley mitotic chromatids. Chromatin from adjacent turns of the helix intermingles due to the stochastic positioning of chromatin loops inside the chromonema. Helical turn size varies along chromosome length, correlating with chromatin density. Constraints on the observable dimensions of sister chromatid exchanges further supports the helical chromonema model.

Evolution of immune genes is associated with the Black Death.

Infectious diseases are among the strongest selective pressures driving human evolution1,2. This includes the single greatest mortality event in recorded history, the first outbreak of the second pandemic of plague, commonly called the Black Death, which was caused by the bacterium Yersinia pestis3. This pandemic devastated Afro-Eurasia, killing up to 30-50% of the population4. To identify loci that may have been under selection during the Black Death, we characterized genetic variation around immune-related genes from 206 ancient DNA extracts, stemming from two different European populations before, during and after the Black Death. Immune loci are strongly enriched for highly differentiated sites relative to a set of non-immune loci, suggesting positive selection. We identify 245 variants that are highly differentiated within the London dataset, four of which were replicated in an independent cohort from Denmark, and represent the strongest candidates for positive selection. The selected allele for one of these variants, rs2549794, is associated with the production of a full-length (versus truncated) ERAP2 transcript, variation in cytokine response to Y. pestis and increased ability to control intracellular Y. pestis in macrophages. Finally, we show that protective variants overlap with alleles that are today associated with increased susceptibility to autoimmune diseases, providing empirical evidence for the role played by past pandemics in shaping present-day susceptibility to disease.

Limitation of current probe design for oligo-cross-FISH, exemplified by chromosome evolution studies in duckweeds.

Duckweeds represent a small, free-floating aquatic family (Lemnaceae) of the monocot order Alismatales with the fastest growth rate among flowering plants. They comprise five genera (Spirodela, Landoltia, Lemna, Wolffiella, and Wolffia) varying in genome size and chromosome number. Spirodela polyrhiza had the first sequenced duckweed genome. Cytogenetic maps are available for both species of the genus Spirodela (S. polyrhiza and S. intermedia). However, elucidation of chromosome homeology and evolutionary chromosome rearrangements by cross-FISH using Spirodela BAC probes to species of other duckweed genera has not been successful so far. We investigated the potential of chromosome-specific oligo-FISH probes to address these topics. We designed oligo-FISH probes specific for one S. intermedia and one S. polyrhiza chromosome (Fig. 1a). Our results show that these oligo-probes cross-hybridize with the homeologous regions of the other congeneric species, but are not suitable to uncover chromosomal homeology across duckweeds genera. This is most likely due to too low sequence similarity between the investigated genera and/or too low probe density on the target genomes. Finally, we suggest genus-specific design of oligo-probes to elucidate chromosome evolution across duckweed genera.

Meiotic crossovers characterized by haplotype-specific chromosome painting in maize.

Meiotic crossovers (COs) play a critical role in generating genetic variation and maintaining faithful segregation of homologous chromosomes during meiosis. We develop a haplotype-specific fluorescence in situ hybridization (FISH) technique that allows visualization of COs directly on metaphase chromosomes. Oligonucleotides (oligos) specific to chromosome 10 of maize inbreds B73 and Mo17, respectively, are synthesized and labeled as FISH probes. The parental and recombinant chromosome 10 in B73 x Mo17 F1 hybrids and F2 progenies can be unambiguously identified by haplotype-specific FISH. Analysis of 58 F2 plants reveals lack of COs in the entire proximal half of chromosome 10. However, we detect COs located in regions very close to the centromere in recombinant inbred lines from an intermated B73 x Mo17 population, suggesting effective accumulation of COs in recombination-suppressed chromosomal regions through intermating and the potential to generate favorable allelic combinations of genes residing in these regions.

Genetic resiliency and the Black Death: No apparent loss of mitogenomic diversity due to the Black Death in medieval London and Denmark.

OBJECTIVES: In the 14th century AD, medieval Europe was severely affected by the Great European Famine as well as repeated bouts of disease, including the Black Death, causing major demographic shifts. This high volatility led to increased mobility and migration due to new labor and economic opportunities, as evidenced by documentary and stable isotope data. This study uses ancient DNA (aDNA) isolated from skeletal remains to examine whether evidence for large-scale population movement can be gleaned from the complete mitochondrial genomes of 264 medieval individuals from England (London) and Denmark. MATERIALS AND METHODS: Using a novel library-conserving approach to targeted capture, we recovered 264 full mitochondrial genomes from the petrous portion of the temporal bones and teeth and compared genetic diversity across the medieval period within and between English (London) and Danish populations and with contemporary populations through population pairwise ΦST analysis. RESULTS: We find no evidence of significant differences in genetic diversity spatially or temporally in our dataset, yet there is a high degree of haplotype diversity in our medieval samples with little exact sequence sharing. DISCUSSION: The mitochondrial genomes of both medieval Londoners and medieval Danes suggest high mitochondrial diversity before, during and after the Black Death. While our mitochondrial genomic data lack geographically correlated signals, these data could be the result of high, continual female migration before and after the Black Death or may simply indicate a large female effective population size unaffected by the upheaval of the medieval period. Either scenario suggests a genetic resiliency in areas of northwestern medieval Europe.

Whole-chromosome paints in maize reveal rearrangements, nuclear domains, and chromosomal relationships.

Whole-chromosome painting probes were developed for each of the 10 chromosomes of maize by producing amplifiable libraries of unique sequences of oligonucleotides that can generate labeled probes through transcription reactions. These paints allow identification of individual homologous chromosomes for many applications as demonstrated in somatic root tip metaphase cells, in the pachytene stage of meiosis, and in interphase nuclei. Several chromosomal aberrations were examined as proof of concept for study of various rearrangements using probes that cover the entire chromosome and that label diverse varieties. The relationship of the supernumerary B chromosome and the normal chromosomes was examined with the finding that there is no detectable homology between any of the normal A chromosomes and the B chromosome. Combined with other chromosome-labeling techniques, a complete set of whole-chromosome oligonucleotide paints lays the foundation for future studies of the structure, organization, and evolution of genomes.

The evolutionary history of dogs in the Americas.

Dogs were present in the Americas before the arrival of European colonists, but the origin and fate of these precontact dogs are largely unknown. We sequenced 71 mitochondrial and 7 nuclear genomes from ancient North American and Siberian dogs from time frames spanning ~9000 years. Our analysis indicates that American dogs were not derived from North American wolves. Instead, American dogs form a monophyletic lineage that likely originated in Siberia and dispersed into the Americas alongside people. After the arrival of Europeans, native American dogs almost completely disappeared, leaving a minimal genetic legacy in modern dog populations. The closest detectable extant lineage to precontact American dogs is the canine transmissible venereal tumor, a contagious cancer clone derived from an individual dog that lived up to 8000 years ago.

Comparative Oligo-FISH Mapping: An Efficient and Powerful Methodology To Reveal Karyotypic and Chromosomal Evolution.

Developing the karyotype of a eukaryotic species relies on identification of individual chromosomes, which has been a major challenge for most nonmodel plant and animal species. We developed a novel chromosome identification system by selecting and labeling oligonucleotides (oligos) located in specific regions on every chromosome. We selected a set of 54,672 oligos (45 nt) based on single copy DNA sequences in the potato genome. These oligos generated 26 distinct FISH signals that can be used as a "bar code" or "banding pattern" to uniquely label each of the 12 chromosomes from both diploid and polyploid (4× and 6×) potato species. Remarkably, the same bar code can be used to identify the 12 homeologous chromosomes among distantly related Solanum species, including tomato and eggplant. Accurate karyotypes based on individually identified chromosomes were established in six Solanum species that have diverged for >15 MY. These six species have maintained a similar karyotype; however, modifications to the FISH signal bar code led to the discovery of two reciprocal chromosomal translocations in Solanum etuberosum and S. caripense We also validated these translocations by oligo-based chromosome painting. We demonstrate that the oligo-based FISH techniques are powerful new tools for chromosome identification and karyotyping research, especially for nonmodel plant species.

Diagnostic microarray for 14 water and foodborne pathogens using a flatbed scanner.

Parallel detection approaches are of interest to many researchers interested in identifying multiple water and foodborne pathogens simultaneously. Availability and cost-effectiveness are two key factors determining the usefulness of such approaches for laboratories with limited resources. In this study, we developed and validated a high-density microarray for simultaneous screening of 14 bacterial pathogens using an approach that employs gold labeling with silver enhancement (GLS) protocol. In total, 8887 probes (50-mer) were designed using an in-house database of virulence and marker genes (VMGs), and synthesized in quadruplicate on glass slides using an in-situ synthesis technology. Target VMG amplicons were obtained using multiplex polymerase chain reaction (PCR), labeled with biotin, and hybridized to the microarray. The signals generated after gold deposition and silver enhancement, were quantified using a flatbed scanner having 2-µm resolution. Data analysis indicated that reliable presence/absence calls could be made, if: i) over four probes were used per gene, ii) the signal-to-noise ratio (SNR) cutoff was greater than or equal to two, and iii) the positive fraction (PF), i.e., number of probes with SNR≥2 for a given VMG was greater than 0.75. Hybridization of the array with blind samples resulted in 100% correct calls, and no false positive. Because amplicons were obtained by multiplex PCR, sensitivity of this method is similar to PCR. This assay is an inexpensive and reliable technique for high throughput screening of multiple pathogens.

Eggshell palaeogenomics: Palaeognath evolutionary history revealed through ancient nuclear and mitochondrial DNA from Madagascan elephant bird (Aepyornis sp.) eggshell.

Palaeognaths, the sister group of all other living birds (neognaths), were once considered to be vicariant relics from the breakup of the Gondwanan supercontinent. However, recent molecular studies instead argue for dispersal of volant ancestors across marine barriers. Resolving this debate hinges upon accurately reconstructing their evolutionary relationships and dating their divergences, which often relies on phylogenetic information from extinct relatives and nuclear genomes. Mitogenomes from the extinct elephant birds of Madagascar have helped inform the palaeognath phylogeny; however, nuclear information has remained unavailable. Here, we use ancient DNA (aDNA) extracted from fossil eggshell, together with target enrichment and next-generation sequencing techniques, to reconstruct an additional new mitogenome from Aepyornis sp. with 33.5X coverage. We also recover the first elephant bird nuclear aDNA, represented by 12,500bp of exonic information. While we confirm that elephant birds are sister taxa to the kiwi, our data suggests that, like neognaths, palaeognaths underwent an explosive radiation between 69 and 52Ma-well after the break-up of Gondwana, and more rapidly than previously estimated from mitochondrial data alone. These results further support the idea that ratites primarily diversified immediately following the Cretaceous-Palaeogene mass extinction and convergently evolved flightlessness. Our study reinforces the importance of including information from the nuclear genome of extinct taxa for recovering deep evolutionary relationships. Furthermore, with approximately 3% endogenous aDNA retrieved, avian eggshell can be a valuable substrate for recovering high quality aDNA. We suggest that elephant bird whole genome recovery is ultimately achievable, and will provide future insights into the evolution these birds.

17th Century Variola Virus Reveals the Recent History of Smallpox.

Smallpox holds a unique position in the history of medicine. It was the first disease for which a vaccine was developed and remains the only human disease eradicated by vaccination. Although there have been claims of smallpox in Egypt, India, and China dating back millennia [1-4], the timescale of emergence of the causative agent, variola virus (VARV), and how it evolved in the context of increasingly widespread immunization, have proven controversial [4-9]. In particular, some molecular-clock-based studies have suggested that key events in VARV evolution only occurred during the last two centuries [4-6] and hence in apparent conflict with anecdotal historical reports, although it is difficult to distinguish smallpox from other pustular rashes by description alone. To address these issues, we captured, sequenced, and reconstructed a draft genome of an ancient strain of VARV, sampled from a Lithuanian child mummy dating between 1643 and 1665 and close to the time of several documented European epidemics [1, 2, 10]. When compared to vaccinia virus, this archival strain contained the same pattern of gene degradation as 20th century VARVs, indicating that such loss of gene function had occurred before ca. 1650. Strikingly, the mummy sequence fell basal to all currently sequenced strains of VARV on phylogenetic trees. Molecular-clock analyses revealed a strong clock-like structure and that the timescale of smallpox evolution is more recent than often supposed, with the diversification of major viral lineages only occurring within the 18th and 19th centuries, concomitant with the development of modern vaccination.

The phylogenetic affinities of the extinct glyptodonts.

Among the fossils of hitherto unknown mammals that Darwin collected in South America between 1832 and 1833 during the Beagle expedition were examples of the large, heavily armored herbivores later known as glyptodonts. Ever since, glyptodonts have fascinated evolutionary biologists because of their remarkable skeletal adaptations and seemingly isolated phylogenetic position even within their natural group, the cingulate xenarthrans (armadillos and their allies). In possessing a carapace comprised of fused osteoderms, the glyptodonts were clearly related to other cingulates, but their precise phylogenetic position as suggested by morphology remains unresolved. To provide a molecular perspective on this issue, we designed sequence-capture baits using in silico reconstructed ancestral sequences and successfully assembled the complete mitochondrial genome of Doedicurus sp., one of the largest glyptodonts. Our phylogenetic reconstructions establish that glyptodonts are in fact deeply nested within the armadillo crown-group, representing a distinct subfamily (Glyptodontinae) within family Chlamyphoridae. Molecular dating suggests that glyptodonts diverged no earlier than around 35 million years ago, in good agreement with their fossil record. Our results highlight the derived nature of the glyptodont morphotype, one aspect of which is a spectacular increase in body size until their extinction at the end of the last ice age.

Combined in vitro transcription and reverse transcription to amplify and label complex synthetic oligonucleotide probe libraries.

Oligonucleotide microarrays allow the production of complex custom oligonucleotide libraries for nucleic acid detection-based applications such as fluorescence in situ hybridization (FISH). We have developed a PCR-free method to make single-stranded DNA (ssDNA) fluorescent probes through an intermediate RNA library. A double-stranded oligonucleotide library is amplified by transcription to create an RNA library. Next, dye- or hapten-conjugate primers are used to reverse transcribe the RNA to produce a dye-labeled cDNA library. Finally the RNA is hydrolyzed under alkaline conditions to obtain the single-stranded fluorescent probes library. Starting from unique oligonucleotide library constructs, we present two methods to produce single-stranded probe libraries. The two methods differ in the type of reverse transcription (RT) primer, the incorporation of fluorescent dye, and the purification of fluorescent probes. The first method employs dye-labeled reverse transcription primers to produce multiple differentially single-labeled probe subsets from one microarray library. The fluorescent probes are purified from excess primers by oligonucleotide-bead capture. The second method uses an RNA:DNA chimeric primer and amino-modified nucleotides to produce amino-allyl probes. The excess primers and RNA are hydrolyzed under alkaline conditions, followed by probe purification and labeling with amino-reactive dyes. The fluorescent probes created by the combination of transcription and reverse transcription can be used for FISH and to detect any RNA and DNA targets via hybridization.

Methods for the preparation of large quantities of complex single-stranded oligonucleotide libraries.

Custom-defined oligonucleotide collections have a broad range of applications in fields of synthetic biology, targeted sequencing, and cytogenetics. Also, they are used to encode information for technologies like RNA interference, protein engineering and DNA-encoded libraries. High-throughput parallel DNA synthesis technologies developed for the manufacture of DNA microarrays can produce libraries of large numbers of different oligonucleotides, but in very limited amounts. Here, we compare three approaches to prepare large quantities of single-stranded oligonucleotide libraries derived from microarray synthesized collections. The first approach, alkaline melting of double-stranded PCR amplified libraries with a biotinylated strand captured on streptavidin coated magnetic beads results in little or no non-biotinylated ssDNA. The second method wherein the phosphorylated strand of PCR amplified libraries is nucleolyticaly hydrolyzed is recommended when small amounts of libraries are needed. The third method combining in vitro transcription of PCR amplified libraries to reverse transcription of the RNA product into single-stranded cDNA is our recommended method to produce large amounts of oligonucleotide libraries. Finally, we propose a method to remove any primer binding sequences introduced during library amplification.

Ancient whole genome enrichment using baits built from modern DNA.

We report metrics from complete genome capture of nuclear DNA from extinct mammoths using biotinylated RNAs transcribed from an Asian elephant DNA extract. Enrichment of the nuclear genome ranged from 1.06- to 18.65-fold, to an apparent maximum threshold of ∼80% on-target. This projects an order of magnitude less costly complete genome sequencing from long-dead organisms, even when a reference genome is unavailable for bait design.

Yersinia pestis and the plague of Justinian 541-543 AD: a genomic analysis.

BACKGROUND: Yersinia pestis has caused at least three human plague pandemics. The second (Black Death, 14-17th centuries) and third (19-20th centuries) have been genetically characterised, but there is only a limited understanding of the first pandemic, the Plague of Justinian (6-8th centuries). To address this gap, we sequenced and analysed draft genomes of Y pestis obtained from two individuals who died in the first pandemic. METHODS: Teeth were removed from two individuals (known as A120 and A76) from the early medieval Aschheim-Bajuwarenring cemetery (Aschheim, Bavaria, Germany). We isolated DNA from the teeth using a modified phenol-chloroform method. We screened DNA extracts for the presence of the Y pestis-specific pla gene on the pPCP1 plasmid using primers and standards from an established assay, enriched the DNA, and then sequenced it. We reconstructed draft genomes of the infectious Y pestis strains, compared them with a database of genomes from 131 Y pestis strains from the second and third pandemics, and constructed a maximum likelihood phylogenetic tree. FINDINGS: Radiocarbon dating of both individuals (A120 to 533 AD [plus or minus 98 years]; A76 to 504 AD [plus or minus 61 years]) places them in the timeframe of the first pandemic. Our phylogeny contains a novel branch (100% bootstrap at all relevant nodes) leading to the two Justinian samples. This branch has no known contemporary representatives, and thus is either extinct or unsampled in wild rodent reservoirs. The Justinian branch is interleaved between two extant groups, 0.ANT1 and 0.ANT2, and is distant from strains associated with the second and third pandemics. INTERPRETATION: We conclude that the Y pestis lineages that caused the Plague of Justinian and the Black Death 800 years later were independent emergences from rodents into human beings. These results show that rodent species worldwide represent important reservoirs for the repeated emergence of diverse lineages of Y pestis into human populations. FUNDING: McMaster University, Northern Arizona University, Social Sciences and Humanities Research Council of Canada, Canada Research Chairs Program, US Department of Homeland Security, US National Institutes of Health, Australian National Health and Medical Research Council.

Quantitative PCR as a predictor of aligned ancient DNA read counts following targeted enrichment.

Targeted DNA enrichment through hybridization capture (EHC) is rapidly replacing PCR as the method of choice for enrichment prior to genomic resequencing. This is especially true in the case of ancient DNA (aDNA) from long-dead organisms, where targets tend to be highly fragmented and outnumbered by contaminant DNA. However, the behavior of EHC using aDNA has been quite variable, making success difficult to predict and preventing efficient sample equilibration during multiplexed sequencing runs. Here, we evaluate whether quantitative PCR (qPCR) measurements of aDNA samples correlate with on-target read counts before and after EHC. Our data indicate that not only do simple target qPCRs correlate strongly with high-throughput sequencing (HTS) data but that certain sample characteristics, such as overall target abundance as well as experimental parameters (e.g., bait concentration and secondary structure propensity), consistently influenced enrichment of our diverse set of aDNA samples. Taken together, our results should help guide experimental design, screening strategies, and multiplexed sample equilibration, increasing yield and reducing the expected and actual cost of aDNA EHC high-throughput sequencing projects in the future.

Detection and characterization of human pathogenic viruses circulating in community wastewater using multi target microarrays and polymerase chain reaction.

Sewage pollution remains the most significant source of human waterborne pathogens. This study describes the detection and characterization of human enteric viruses in community wastewaters using cell culture coupled with multiple target microarrays (with a total of 780 unique probes targeting 27 different groups of both DNA and RNA viruses) and polymerase chain reaction (PCR) assays. Over a 13-month sampling period, RNA viruses (astroviruses and enteroviruses) were more frequently detected compared to DNA viruses (adenoviruses, particularly type 41 and BK polyomavirus). Overall, many more viruses were shed during the winter months (December-February) compared to the summer months. Exploration of the multiple types of enteric viruses particularly in winter months identified much more significant prevalence of key viral pathogens associated with sewage pollution of the water environment than previously realized and seasonal disinfection used in some parts of the world may lead to a seeding of ambient waters. Molecular characterization of pathogenic viruses in community wastewater will improve the understanding of the potential risk of waterborne disease transmission of viral pathogens.

From design to screening: a new antimicrobial peptide discovery pipeline.

Antimicrobial peptides (AMPs) belong to a class of natural microbicidal molecules that have been receiving great attention for their lower propensity for inducing drug resistance, hence, their potential as alternative drugs to conventional antibiotics. By generating AMP libraries, one can study a large number of candidates for their activities simultaneously in a timely manner. Here, we describe a novel methodology where in silico designed AMP-encoding oligonucleotide libraries are cloned and expressed in a cellular host for rapid screening of active molecules. The combination of parallel oligonucleotide synthesis with microbial expression systems not only offers complete flexibility for sequence design but also allows for economical construction of very large peptide libraries. An application of this approach to discovery of novel AMPs has been demonstrated by constructing and screening a custom library of twelve thousand plantaricin-423 mutants in Escherichia coli. Analysis of selected clones by both Sanger-sequencing and 454 high-throughput sequencing produced a significant amount of data for positionally important residues of plantaricin-423 responsible for antimicrobial activity and, moreover, resulted in identification of many novel variants with enhanced specific activities against Listeria innocua. This approach allows for generation of fully tailored peptide collections in a very cost effective way and will have countless applications from discovery of novel AMPs to gaining fundamental understanding of their biological function and characteristics.